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Abstract 


The  fractakdi  mention  of  a  curve  is  a  measure  of  its  geometric  complexity  and 
can  be  any  nbn-integer  value  between  1  and  2  depending  upon  the  curve's  level 
of  complexity.  ^This  paper  discusses  an  algorithm,  which  simulates  walking  a 
pair  of  dividers  along  a  curve,  used  to  calculate  the  fractal  dimensions  of 
curves.  It  also  discusses  the  choice  of  chord  length  and  the  number  of  solution 
steps  used  in  computing  fracticality.  Results  demonstrate  the  algorithm  to  be 
stable  and  that  a  curve’s  fractal  dimension  can  be  closely  approximated. 
Potential  applications  for  this  technique  Include  a  new  means  for  curvilinear 
data  compression,  description  of  pianimetric  feature  boundary  texture  for 
improved  realism  in  scene  generation  and  possible  two-dimensional  extension 
for  description  of  surface  feature 


INTRODUCTION 


The  problem  of  describing  the  forms  of  curves  has  vexed  researchers  over  the 
years.  For  example,  a  coastline  is  neither  straight,  nor  circular,  nor  elliptic  and 
therefore  Euclidean  lines  cannot  adquataly  describe  most  real  world  linear 
features.  Imagine  attempting  to  describe  the  boundaries  of  clouds  or  outlines 
of  complicated  coastlines  In  terms  of  classical  geometry.  An  intriguing  concept 
proposed  by  Mandelbrot  (1967, 1977)  is  to  use  fractals  to  fill  the  void  caused  by 
the  absence  of  suitable  geometric  representations.  A  fractal  characterizes 
curves  and  surfaces  in  terms  of  their  complexity  by  treating  dimension  as  a 
continuum.  Normally,  dimension  Is  an  briefer  number  (1  for  curves,  2  for  areas, 
and  3  for  volumes))  however,  fractal  dimensions  may  vary  anywhere  between  1 
and  2  for  a  curve  and  2  and  3  for  a  surface  depending  upon  the  Irregularity  of 
the  form.  Although  individual  fractals  have  been  around  since  the  1909s  , 

first  to  reconnlae  their  applications  outride  of 


This  paper  dUcmw  an  algorithm,  writttn  In  an  interactive  setting,  dwipwd  to 
measuro  tha  fractkality  of  a  cun*  and  additions  to  theory,  it  also  presents 
results  from  examining  several  cartographic  curves. 

DEFINITION  OF  FRACTALS  AND  SELF-SIMILARITY 

In  Euclidean  geometry  every  curve  has  a  dimension  of  1  and  every  plane  has  a 
dimension  of  2.  This  is  generally  referred  to  as  the  topological  dtotenstan  (DtL 
These  dimensions  remain  constant  no  matter  how  complex  or  irregidar  a  curve 
or  plane  may  be.  For  example,  the  west  coast  of  Greet  Britain  contains  many 
irregularities,  but  the  topological  dimension  remains  I. 

In  the  fractal  domain  a  curve's  dimension  may  be  between  1  and  2  according  to 
its  complexity.  The  more  contorted  a  straight  line  becomes,  the  higher  its 
fractal  dimension.  Similarly,  a  plane's  dimension  may  be  a  non-integer  value 
between  2  and  3.  The  fractal  dimension  for  any  curve  or  surface  Is  denoted  by 
(D)  and  within  this  framework:  D  >Dt.  Mandelbrot  (1977)  proposes  the 
following  definition  for  a  fractal:  MA  "Fractal  will  be  defined  as  a  set  for  which 
the  Hausdorff-Besicovitch  dimension  strictly  exceeds  the  topological  dimen¬ 
sion."  The  precise  definition  of  the  Hausdorff-Besicovitch  dimension  can  be 
found  in  Besicovitch  and  Ursell  (1937). 

Suppose  one  measures  the  length  of  a  curve  using  dividers  set  to  a  prescribed 
step  size  (n)  and  counts  the  number  of  steps  (N)  it  takes  to  approximate  the 
length.  If  a  curve  is  a  fractal,  then  its  length  approaches  infinity  as  (n)  goes  to 
zero.  With  this  nomenclature,  Mark  (1979;  states  the  Hausdorff-Beeocoviteh 
dimension  of  a  curve  is  the  absolute  value  of  the  slope  (D)  of  a  power  function 
relating  N  (n)  to  n  where  n  is  a  step  size  and  the  number  of  steps  along  the 
curve,  N  (n)  is  counted.  In  summary,  D  is  a  measure  of  the  complexity  for  any 
curve  or  surface  and  all  fractal  sets  are  nowhere  differentiable. 

Central  to  the  concept  of  fractals  is  the  notion  of  self -similarity.  Self- 
similarity  means  that  for  any  curve  or  surface  a  portion  of  the  curve  or  surface 
can  be  considered  a  reduced  Image  of  the  whole.  However,  seldom  In  nature 
(crystals  are  one  exception)  does  self -similarity  occur  and  therefore  a  statisti¬ 
cal  form  of  self-similarity  Is  often  encountered.  In  other  words,  if  a  curve  or 
surface  is  examined  at  any  scale  it  will  resemble  the  whole  In  a  statistical 
sense;  therfore,  D  will  remain  constant.  Brownian  motion  Is  an  excellent 
example  of  statistical  self -similarity.  Because  of  this  principle,  a  curve  can  be 
decomposed  into  N=r  nonoverlapping  parts  and  each  subsegment  has  a  length  of 
1/rsl/N.  Similarly,  a  unit  square  can.be  efivided  into  N*r*  squares,  where  the 
similarity  ratio  is  KN)  *  1/r  *  1/N  .  In  either  case  the  following  equation 
applies: 

D«  log  N/log  (1/r)  (1) 

and  could  be  called  the  shape’s  similarity  dimension.  D  can  also  be  expressed 
as: 


D*log  (N/N0Vk>g  <VA)  (2) 

where  A*  and  X  are  two  sampling  intervals  and  N  and  N  are  the  number  of 
such  intervals  contained.  If  a  curve  resembles  a straight  lint  then  when  the 
sampling  interval  Is  halved,  N  doubles  and  the  proportion  equals  1.  The 
majority  of  cartographic  curves  are  not  straight  lines  and  therefore  N  will  more 
than  double  causing  D  to  be  greater  then  1.  The  principle  of  seii-MffiilarUy  is 
dbmiseed  by  GooddtUd  (ltio),  Hakanoon  (1971),  and  SchtMsggtr  (1970). 
Hakanson,  for  example,  points  out  the  abeurdlty  of  postuSatk^  tha  validity  of 


Richardson  (1961)  found  the  west  coast  of  Britian  to  have  a  constant  D  of  1*25 
over  sampling  intervals  between  10  and  I000km„  he  found  the  east  coast  to 
vary  between  1,19  and  1J1  for  a  similar  sampling  interval.  This  suggests  that 
whatever  created  the  irregularities  on  the  coastline  acted  at  specific  scales. 
Goodchild  states  that  since  seif -similarity  is  only  one  aspect  of  the  fractal 
approach,  it  would  be  unwise  to  reject  the  entire  concept, 

DEVELOPMENT  OF  THE  FRACTAL  CURVE  ALGORITHM 
AND  EXTENSION  OF  THEORY 

The  following  original  algorithm  is  based  on  the  earlier  empirical  work 
performed  by  Richardson  (1961)  and  later  extended  by  Mandelbrot  (1967), 
Richardson  measured  the  lengths  of  several  frontiers  by  manually  walking  a  pair 
of  dividers  along  the  outline  so  as  to  count  the  number  of  steps.  The  opening  of 
the  dividers  (n)  was  fixed  in  advance  and  a  fractional  side  was  estimated  at  the 
end  of  the  walk.  The  main  purpose  in  this  section  of  Richardson's  research  was 
to  study  the  broad  variation  of  2h  with  n. 

Richardson  produced  a  diagram  in  which  he  plotted  log  total  length  gainst  log 
step  size  and  it  is  shown  in  Figure  i.  Mandelbrot  (1967)  discovered  a 
relationship  between  the  slope  (0)  of  the  lines  and  fractal  dimension  (D).  To 
Richardson  the  slope  had  no  theoritical  meaning,  but  to  Mandelbrot  It  could  be 
used  as  an  estimate  of  1-D,  which  leads  tot 


D*  1-0  “  (3) 

The  algorithm  simulates  walking  a  pair  of  dividers  along  a  curve  and  counts  the 
number  of  steps.  In  cases  where  more  than  one  intersection  occurs,  the 
intersection  which  comes  first  in  order  forward  along  the  curve  is  selected.  To 
be  more  accurate,  step  size  (prescribed  opening  of  the  dividers)  is  called  chord 
length  (cl)  and  the  number  of  steps  is  called  the  number  of  chord  lengths. 

In  order  to  begin  walking  the  dividers  along  the  curve,  the  dividers  must  be  set 
to  some  opening.  The  curves  used  in  this  research  are  not  infinitely  subdivided 
fractal  curves  so  that  selection  of  the  Initial  chord  length  must  be  based  on 
some  attribute  of  the  curve.  For  a  very  contorted  curve  it  would  be 
meaningless  to  choose  a  chord  length  many  times  smaller  than  the  smallest  line 
segment.  If  an  extremely  small  chord  length  is  selected,  an  attempt  to  examine 
the  fractal  character  of  a  curve  would  extend  beyond  the  primitive  subeiements 
used  to  represent  the  geometry  of  the  resulting  form.  In  other  words,  beyond 
this  lower  limit  of  primitive  subelements,  the  curve's  fractal  dimension  behaves 
as  if  it  is  a  straight  line.  A  suggested  initial  chord  length  is  determined  by 
calculating  the  cfistance  between  each  two  consecutive  points  on  the  curve  and 
taking  1/2  the  average  distance.  The  average  <fistance  is  divided  by  2  because 
the  sampling  theorem  states  one  should  sample  at  1/2  the  wavelength  so  that  no 
significant  variation  escapes.  This  presents  an  approximate  lower  limit  as  to 
the  selection  of  the  initial  chord  length.  Although  the  accuracy  of  this  method 
is  dependent  on  the  manner  in  which  the  curve  Is  digitized,  the  form  of  the 
curve  often  dictates  this  manner. 

After  the  initial  chord  length  is  determined,  the  algorithm  computes  the 
distance  between  the  first  two  points  on  the  curve  using  the  standard  distance 
formula.  If  the  dfoance  is  greater  than  chord  length  (d)»  a  new  point  is 
interpolated  between  points  I  and  2  using  the  following  interpolation  equations: 

DP  *  (d-DZSTl)  /  (D1ST1-DISTA)  (6) 

XNEW  *  X!  ♦  DP*  (XrXj)  (3) 

YNEW  *  Y|  +DP*  (Y 2*Y  j)  (6) 


****** 


Jl 


where  DP  *  dUtance  proportion 

DISTi  *  distance  between  the  present  point  and  the  first  forward 
point  on  the  curve 

DISTA  a  distance  between  the  present  point  and  the  second  forward 

point  on  the  curve 

XNCV  *  new  X-  coordinate 

YNEW  «  new  Y-  coordinate 

X9Y  3  X  and  Y  coordinates  of  point  1  and  2* 

Figure  2  demonstrates  how  a  point  is  interpolated  on  a  straight  line  segment* 

If  the  distance  is  less  than  the  chord  length,  the  dbtance  between  points  1  and  3 
(DISTC)  is  computed*  If  D1STC  is  greater  than  the  chord  length,  it  is  known 
that  the  chord  length  segment  intersects  between  points  2  and  3  and  that  the 
distance  between  these  points  is  determined  (DISTB);  See  Figure  3a*  The  point 
of  intersection  is  computed  using  trigonometric  functions.  An  angle  C  is  deter¬ 
mined  using  the  law  of  cosines* 

C.  cos'1  DISTB2  ♦  DISTA2  -  DISTC2  (7) 

FtWSTB  •  DISTA 

Since  angle  C  is  known,  an  angle  A,  which  is  the  ai«le  the  chord  length 
intersects  between  points  2  and  3,  can  be  computed. 

A  =  SIN'1  «DISTA*sin  C)/d)  («) 

Now  that  two  angles  are  known,  angle  B  is  easily  computed.  Because  angles  A 
and  B  are  known,  a  side  (DISTBO  can  be  calculated;  see  Figure  3b. 

DISTB1  =  DISTA*sin  B/sin  B  (9) 

DISTB’  provides  the  distance,  from  point  2,  in  which  the  chord  length's 
intersection  is  located  on  the  segment  between  points  2  and  3.  A  distance 
proportion  is  calculated  usings 

DP  r  DISTB'/DISTB  (10) 

Since  the  distance  proportion  and  the  X,Y  coordinates  for  points  2  and  3  are 
known,  the  equations  used  to  interpolate  for  a  straight  line  segment  can  be  used 
to  determine  the  new  coordinates;  see  Figure  3c.  After  the  new  point  is 
located,  this  new  point  becomes  point  1  and  the  next  two  forward  points  on  the 
curve  become  points  2  and  3.  Each  time  a  chord  length's  intersection  is 
determined,  1  is  added  to  the  number  of  chord  lengths. 

In  the  case  where  DISTC  Is  less  than  the  chord  length,  the  third  point  is 
incremented  by  1  (fourth  point)  and  the  dbtance  again  checked.  This  continues 
until  the  distance  is  greater  than  the  chord  length  or  the  end  of  the  curve  is 
encountered;  see  Figure  h  When  the  distance  does  become  greater  than  the 
chord  length,  the  chord  length's  point  of  intersection  is  determined  by  using  the 
same  trigonometric  equations  as  discussed  above.  The  only  difference  is  the 
sides  of  the  tthngles  may  be  loiter.  At  the  end  of  the  curve,  if  the  chord 
length  is  greater  than  DISTA,  the  portion  of  the  remaining  chord  length  is  added 
to  the  number  of  chord  lengths. 

After  the  dfoiders  are  walked  along  the  curve  with  the  initial  chord  iewth,  the 
dividers  are  opened  to  another  distance.  This  distance  Is  a  geometric  adding  of 
the  first  chord  length.  For  example,  If  the  Initial  chord  la*th  Is  2,  then  the 
subsewient  chord  lengths  would  be  t,  S,  16,  32,  66,  and  so  on.  This  eliminates 
biasing  when  using  linear  regression  because  on  a  logarithmic  scale,  geometric 
adding  provides  equal  spacing  between  the  chord  length*. 


FHniTTT'T 


•  i.ai« 


M»w«  S.  Kodiak  Ialaad-163)  polata. 


The  number  of  solution  steps  or  the  number  of  times  the  dividers,  with  different 
chord  lengths,  ere  walked  along  the  curve  is  limited  by  the  number  of  chord 
lengths  used  to  estimate  length.  The  minimum  number  of  chord  lengths  used  to 
approximate  length  is  3.  This  is  chosen  to  provide  consistency  among  results  as 
opposed  to  using  a  variable  limit,  but  is  subject  to  change  pending  additional 
research. 

After  each  time  the  dividers  are  walked  along  the  curve,  the  number  of  chord 
lengths  and  the  corresponding  chord  lengths  are  saved.  These  are  used  in  the 
linear  regression  where  log  line  length  (number  of  chord  lengths  *chord  length) 
is  regressed  against  log  chord  length.  A  curve’s  fractal  dimension  is  determined 
by  using  equation  3. 

To  provide  an  indication  of  the  proportion  of  variance  in  the  response  variable 
explained  by  the  describing  variable,  r*  is  computed.  This  value  plays  ap 
important  role  in  determining  the  optimum  number  of  solution  steps.  A  low  r  , 
for  example  when  the  number  of  solution  steps  equals  12, jneans  the  initial 
chord  length  falls  below  the  primitive  subelements.  A  low  r  is  determined  by 
decreasing  the  number  of  solution  steps  and  comparing  the  two  values.  The 
desirable  number  of  solution  steps  is  indicated  when  r  reaches  its  maximum 
without  the  number  of  steps  failing  below  3.  A  linear  regression  with  less  than 
3  points  opens  up  some  criticisms  as  to  the  validity  of  results  and  it  should  be 
emphasized  the  linear  regression  is  used  as  a  parameter  estimate  and  not  for 
statistical  inferences. 

EXAMPLES  AND  RESULTS 

Of  the  five  land-frontiers  Richardson  measured,  low  have  been  point  digitized 
and  used  in  this  study.  They  are:  coast  of  the  Australian  mainland;  coast  of 
South  Africa,  starting  from  Swakopmund  to  Cape  Sta.  Lucia;  frontier  between 
Spain  and  Portugal,  moving  south  to  north;  and  the  west  coast  of  Great  Britain, 
from  Land's  End  to  Duncansby  Head.  Table  1  shows  D  as  the  result  of 
Richardson's  measurements  and  the  new  D  suggested  by  this  research.  The 
expected  discrepancy  is  the  result  of  the  digitization  process  because  digitiza¬ 
tion  allows  the  capture  of  minute  detail  in  a  curve,  and  since  these  curves  were 
digitized  at  a  larger  scale,  a  higher  D  is  anticipated. 


Curve 

Slone  (B) 

New  D 

West  Coast  of  Great  Britain 

-.23 

1.23 

1.2671 

The  Coast  of  the  Australian 
Mainland 

-.13 

1.13 

1.1490 

Cost  of  South  Africa 

-.02 

1.02 

1.0336 

Land-frontier  between  Spain 
and  Portugal 

-.14 

1.14 

1.1014 

Table  1.  Result  from  Richardson's  (i960  research,  corresponding  fractal 
dimension  and  the  new  suggested  fractal  dimension. 


For  this  paper,  Kodiak  bland  is  taed  to  demonstrate  how  the  fractal  curve 
algorithm  operates.  The  curve  w&  digitized  in  trace  mode  with  delta  minimum 
and  delta  maximum  variations  at  .002  and  .03  respectively.  The  outline 
contains  1633  points  and  is  in  Figure  3.  The  results  from  calculating  D  are  in 
Table  2  where  the  different  initial  chord  lengths  are  selected  to  show  the 
possible  variations  in  D  over  a  number  of  sampling  intervals.  The  results  show 
D  to  vary  from  1.1S36  to  1.3719.  These  variations  in  D  reflect  a  lack  of  seif- 
similarity  in  the  curve. 

The  selection  of  an  extremely  short  initial  chord  length  of  .01S33  represents 
examining  the  curve  below  its  primitive  subelements  and  biases  D  toward  a 


llllill 
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straight  line.  The  corresponding  scatterplot  in  Figure  6  displays  a  curvature  of 
the  data  points  resulting  in  a  lower  R-SQ.  It  is  this  type  of  curvature, 
resembling  the  shape  of  a  rainbow,  that  indicates  the  shortness  of  the  chord 
length.  The  chord  length  of  .01833  is  the  average  segment  length  and  is 
calculated  by  computing  the  distance  between  each  two  consecutive  points  on 
the  curve,  summing  the  distances,  and  dividing  by  the  number  of  segments. 
The  suggested  D  to  represent  the  curve  is  1.3103.  The  most  appropriate  D 
value  is  determined  from  the  minimum  amount  of  curvature  present  in  the 
scatterplot  resulting  in  the  relatively  high  R-SQ  value.  The  suggested  initial 
chord  length  of  .03666  is  still  too  small,  indicated  by  the  low  R-SQ  value,  but 
represents  a  starting  point  at  which  to  determine  D. 

A  thinned  version  of  Kodiak  Island  is  in  Figure  7  and  contains  1000  points.  The 
elimination  of  633  points  is  accomplished  with  a  program  which  deletes 
excessive  points  within  a  certain  chord  length.  Table  3  indicates  D  varying 
between  1.214  and  1.3639.  The  comparison  between  the  same  chord  length  of 
.05894  for  the  original  and  thinned  islands  displays  how  stable  the  algorithm  is 
to  measure  D.  This  initial  chord  length  produced  a  D  of  1.3105  (1633  points) 
and  1.2949  (1000  points)  giving  a  1.19%  difference.  The  D  for  the  1000  point 
island  is  expected  to  be  slightly  lower  because  any  data  thinning  process 
normally  removes  some  complexity  from  the  feature.  The  proposed  D  for  the 
thinned  island  is  approximately  1.2949  and  the  scatterplot  is  in  Figure  8. 

SUMMARY  AND  CONCLUSIONS 

The  results  based  on  the  previous  empirical  curves  point  out  the  importance  of 
selecting  the  appropriate  initial  chord  length.  A  chord  length  which  is  too  short 
is  easily  detected  by  either  examining  the  amount  of  curvature  present  in  the 
scatterplot  or  the  low  R-SQ  value.  Normally,  the  suggested  initial  chord  length 
falls  within  this  category,  but  it  must  be  emphasized  that  this  chord  length  is 
merely  a  beginning  point.  The  ideal  initial  chord  length,  which  produces  the 
most  appropriate  D,  is  selected  by  observing  the  behavior  of  the  scatterplo ts, 
R-SQ  values,  and  the  solution  steps.  This  research,  like  Richardson’s  work, 
indicates  that  from  5  to  8  solution  steps  are  sufficient  to  determine  the  slope  of 
the  regression  line  and  thus  fracticality. 

The  results  also  indicate  the  fractal  curve  algorithm  to  be  stable,  and  that  it  is 
able  to  closely  approximate  D.  The  variations  in  D,  over  a  number  of  sampling 
intervals,  reflect  a  need  to  examine  the  effects  of  seif-similarity,  or  lack  of  it 
on  a  curve’s  fracticality.  Finally,  this  research  brings  into  focus  the  strong 
problem  solving  capabilities,  at  the  hands  of  cartographers,  through  the  use  of 
interactive  computer  graphics. 
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